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We investigate some of the experimental, observational and theoretical consequences of hypo- 
thetical stable black holes in the mass range between the electro-weak scale and the Planck mass, 
2.4x10^^ TeV. For the purpose of calculations we use Lovelock black holes in odd dimensions. If 
such black holes exist they contribute to dark matter. We show that the passage of the black holes 
through matter and the collision of black holes have a well defined experimental signature. Depend- 
ing on their cross section and energy they also accumulate in stars and influence their development. 
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INTRODUCTION 

Gravity theories in space-time dimensions greater than four have been studied extensively ever since the seminal 
papers of Kaluza and Klein on unification of gravity and gauge forces were published in the early twentieth century. 
Modern attempts at unification of forces at the quantum level, supergravity and super string theory are formulated 
in eleven and ten dimensions respectively. Low scale gravity models proposed by Arkani Hamed, Dimopolous and 
Dvali [1] or Randall and Sundrum [2] [2] to explain the hierarchy problem are also formulated in dimensions greater 
than four. In all these models bulk extra dimensions have to be compactified or made small by warping to explain 
why the observed universe has only four space-time dimensions. 

These extra dimensional gravity models admit black hole or black brane solutions. Masses of these black objects 
could range from the fundamental Planck scale of the theory up to astrophysical mass scales of millions of solar 
masses. Microscopic black holes in these models would have sizes smaller than the scale of compactification while 
large astrophysical black objects would more likely be black strings or branes. Planck scales in these models could 
vary from standard Mp = 2.4 x 10^^ TeV down to 1 TeV, in low scale gravity models. In particular, low scale gravity 
models admit microscopic black holes of TeV mass. Such objects could be produced in particle accelerators like the 
LHC [4 or in collisions of high energy cosmic rays with matter [S]. Most of the studies investigating black hole 
production at the LHC assume that the produced black holes are of the standard Schwarzschild-Tangherlini type and 
will decay almost instantaneously into elementary particles via Hawking radiation [6] . An exception is the recent work 
of Giddings and Mangano 7 , which examines the consequences of the very conservative assumption that Hawking 
radiation does not exist at all and conclude that TeV scale black holes produced at the LHC would pose no risk to 
the Earth over its lifetime. 

In this paper we will investigate the fate of the black hole relics in a model, which incorporates standard Hawking 
radiation but in which black hole masses have a lower bound. Upon approaching the mass bound the Hawking 
temperature of these black holes goes to zero and they become stable in vacuum. Their cross section for collisions 
with elementary particles and each other would, however, stay finite, leading to observable consequences when they 
come into contact with matter or collide with each other inside stars or in the galaxy. 

If there are stable microscopic primordial black holes then one has to confront the possibility that they could have 
a relic density large enough to over close the universe. For this to become a serious issue the temperature of the 
universe at some epoch has to be comparable to the minimum stable mass of the black holes so that they could be 
produced in particle collisions. For the current discussion we assume that the reheating temperature of the universe 
after inflation is always much less than the minimum mass of the black holes. We plan to discuss the issue of the relic 
density of stable Lovelock black hole remnants in another publication. 

According to our hypothesis all primordial black holes, or black holes produced by cosmic rays, that would otherwise 
have decayed by the present epoch should end up in (or very close to) the minimal mass end state. Such black holes 
contribute to or may even dominate dark matter. Their number can only decrease in collisions with each other or by 
absorption by macroscopic black holes. We will discuss the passage of stable black holes through matter and point 
out that their interaction with matter would have a characteristic signature, detectable in underground experiments. 

Lovelock gravity theories admit black hole solutions with a limiting minimum mass and vanishing Hawking tem- 
perature. There is no reason known to us why in dimensions higher than four the contribution of Lovelock terms [9], 
formed from higher powers of the curvature tensor, to the Einstein action should be suppressed. Now in odd dimen- 
sions the ADM masses of black holes are limited below by the coupling strength of Lovelock terms. Furthermore, 
when such black holes evaporate by Hawking radiation their Hawking temperature approaches zero when their ADM 
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mass approaches the lower bound. The prime example of this is black holes in 5-d Einstein Gauss-Bonnet gravity. 
Therefore, such a scenario is not completely in the realm of phantasy. 

Since the world has 4 non-compact dimensions on macroscopic scales, we need to compactify all but four of the 
odd number of dimensions in which the stable minimum mass black holes live. If the size of a black hole is smaller 
than the compactification scale, then the caged black hole should still have a stable minimum mass. A systematic 
expansion in the inverse of the compactification scale for black holes in Einstein gravity has been studied in the past 
[25] Lovelock black holes can be embedded into a higher dimensional space, turning the black holes into black 
branes in those extra dimensions. While we return to the discussion of these questions in the last section of this paper 
we will study these and similar problems more extensively in a future publication. 

A simpler alternative, of trying to compactify all dimensions beyond 4 at the same scale in Gauss-Bonnet gravity 
has recently been attempted |10j . Though a minimal mass exists in such a scenario, but the Hawking temperature 
does not vanish and consequently black branes are not stable at the minimal mass. 

The consequences of Gauss-Bonnet gravity on black hole phenomenology have already been considered by 
Rizzo [11] [12] concentrating on the signature of TeV-scale stable black holes produced at LHC, and by Alexeyev, 
Barrau, Boudoul, Khovanskaya, and Sahzin [13] in the framework of an ADD [T] or Randall-Sundrum II [5] model 
with limiting masses of macroscopic range (1 g- 10^* g). We would like to concentrate on models in which the minimal 
mass of black holes is comparable or below the 4 dimensional Planck mass. Such a mass limit would be natural for 
example in the framework of a Randall-Sundrum [2] scenario, ADD theory in 9 or 11 dimensions, or Kaluza-Klein 
models in 11 dimensions and some string-inspired models of compactification. 

Though Schwarzschild-Lovelock black holes in odd dimensions are at the back of our minds, our investigations are 
fairly independent of specific models. They require the knowledge of a few basic parameters, namely the coupling 
strength of the Lovelock terms, to describe the near-minimal mass behavior of these solutions. We assume that all 
non-topological Lovelock terms are present in the Lagrangian and all of their couplings are in the scale of the Planck 
mass. It is worth mentioning that higher order corrections to heterotic string theories contribute a Gauss-Bonnet 
term (second order Lovelock term) to the classical action [M], jl5j . 

When all the Lovelock coefhcients are of 0(1), black holes or black branes with M >> Mp are well described by 
the Schwarzschild-Tangherlini metric, their behavior is quite standard and well described by Einstein gravity. In this 
paper we will concentrate on black holes with masses close to the minimum mass of the order of Mp, when Lovelock 
terms play an important role leading to significant phenomenological consequences. 

In the current paper we have not addressed questions pertaining to detection of stable microscopic black hole 
remnants by direct dark matter detection experiments [16] . |17j . On one hand, if the minimum mass of these black 
holes is close to 1 TeV, then we will show that their abundance must be extremely low compared to other components 
of dark matter to avoid a contradiction with proton decay experiments. On the other hand, if their mass is much 
higher it could be difficult to detect them in direct dark matter detection experiments. 

Finally, to close this section we list the issues we intend to address in the rest of this paper, (i) Stable minimal mass 
black holes aggregate in clumps under the influence of gravity. Then they will form part (or all) of dark matter. This 
effect was considered for macroscopic stable black holes by Alexeyev et. al. [T3], and for even larger, astronomical 
size black holes by Frampton [18j . as well. These black holes may accumulate inside stars depending on the parton 
cross section, their velocity distribution, and the gravity at the surface of the star. We will discuss the rate of heat 
production by the collision of accumulated black holes and compare it to that of "standard" WIMPs [19]. (ii) Black 
holes, as part of dark matter, would collide with partons while passing through detectors in underground laboratories. 
These collisions would be followed by a characteristic decay that is quite different from those of Schwarzschild black 
holes. Using past proton decay experiments we are able to give a rough lower bound on the minimum mass of stable 
black holes and/or their abundance compared to other dark matter constitutents. We also propose an experimental 
investigation of finding stable black holes in underground laboratories, (iii) Relativistic black holes produced by cosmic 
rays [5] passing through stars or planets may evaporate to minimal size or increase in size by accretion depending on 
the interaction cross section. We estimate what is the lowest energy of incoming cosmic rays that allows black holes 
to increase in size and accrete all available matter in a star. 

LOVELOCK BLACK HOLES IN ODD DIMENSIONS 

As we mentioned in the introduction our investigations are largely independent of the choice of a specific gravity 
model. We only assume that black holes have a minimal mass, at which the Hawking temperature vanishes. Lovelock 
black holes [9] in odd dimensions are prime examples satisfying these constraints. Gravitational theories including 
Lovelock terms are perfectly acceptable, as their field equations contain only second order derivatives and satisfy 
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standard requirements on a classical theory of gravity just as much as Einstein gravity does. The langrangian of the 
nth order Lovelock term, L„, is built from an antisymmetrized combination of the nth power of the curvature tensor. 
In D dimensions a linear combination of Lovelock terms up to n = (D — l)/2, such as 

{D-l)/2 

L = R+ J2 a4^""^)«„L„ (1) 

n=2 

contribute non-trivially to the equations of motion. Note that n = corresponds to the cosmological constant, the 
n = 1 term is Einstein gravity, and Gauss-Bonnet gravity corresponds to n = 2. We extract an appropriate power of 
the Planck mass to make the coupling constants a„ dimensionless. Also note that we use Mp as a unknown quantity, 
dependent on D, to be determined by experiments. We also keep the coupling constants, a„ as free parameters, but 
when we need to make an estimate of an observable quantity we will often consider predictions under the assumption 
o^n = 0{l), or, alternatively Mmin ~ Mp, where Mj„in is the minimal mass. We use a definition of Newton's constant 
in terms of the Planck mass of the "Particle Data Book" 121j, 

SttGd = {27t)^'HiI~^ (2) 

With this definition the experimental lower bound on Mp is Mp > 1 TeV [3T]. This bound was calculated assuming 
Schwarzschild-Tangherlini black holes. We expect that the lower bound on the Planck mass when black holes have a 
minimal mass is somewhat lower. The reason for that is the small amount of radiation produced when the black holes 
approach the minimal mass. Since it barely radiates one could only recognize its production by the missing energy 
that could have been taken away by other neutral particles. 

Analytic spherically symmetric black hole solutions have been found for n = 2 [Tn][12][lS] and 3 [23] Lovelock 
gravity. While in even dimensions black hole masses have no lower bound, in odd dimensions the mass is limited by 
the Lovelock coupling a„ , to ensure the existence of a horizon. In particular for D ~ 5 and 7 

(Ai(^)Mp)^-3 >{D^ 'S)la^D-i)/2- (3) 

In ([3]) /i^'^-' is the radius of horizon for the D-dimcnsional Tangherlini black hole [7] 

where 



2(27r 



,D-4 



{D-2)n(D-2)' 



(5) 



and 



-r[p-i)/2]- 



Using the latter equations the lower bounds on the ADM mass are 

M > M^l^37ra2Mp, 



M > Mj^l = yaaMp (7) 

for D = 5 and D = 7, respectively. It is worth noting that the minimal mass is larger than the Planck mass if the 
coupling of the Lovelock terms, a„ = 0(1). This result may be of significance if we consider the stability of minimum 
masses under quantum gravity. It is possible that quantum effects destabilize classically stable black holes. In what 
follows, we assume that even if that were true, their lifetimes are long enough to leave our classical calculations valid. 
The radius of horizon is 



rf = Af^y (Ai(5)Mp)2 - 2^2, 



= Afp^/-6a2 + V0i(7)Mp)4-24a3 + 36ai, (8) 
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in 5 and 7 dimensions, respectively. 

The Hawking temperature in 5 dimensions is 

In D=7 the expression for Th is much more complicated, but also proportional to the radius of horizon, as in D = 5. 
The reason for this is simple. In odd dimensions the metric functions are even functions of the radius and the surface 
gravity is proportional to dgtt/dr. Then both of these theories the radius of horizon and the Hawing temperature 
vanish as ~ \JM — Mmin . This is a general feature of theories we consider in this paper. 

Strictly speaking minimum mass black holes do not exist. Scalars, like the curvature scalar, square of the Ricci 
tensor, and the Kretschmann scalar all diverge at zero radial coordinate, when the radius of the horizon vanishes. 



This state is, however, never reached in the semi-classical theory because, as is shown by (13), it takes infinite time 
for the mass of a black hole produced with M > Mmin to reach Mmin ■ 

While at sufficiently large mass, M > Mp, the Hawking temperature approaches its value in Einstein gravity 
Tff ~ (27rr;j)^^, at M ~ Afmin in Lovelock gravity it can be parametrized by the constants, Mp and q;„ as 

-'-H - 




V mm 

In general we may parametrize the dependence of the Hawking temperature on the minimal mass by the two parameters 
Ms and Mp as 



/ M 

V JWmin 

where Ms,Afmin = 0{Mp). We will use this formula in all subsequent calculations. However, whenever we need 
concrete numbers for calculating a particular phenomenon related to stable black holes we will rely on the 5 dimensional 
formulas, which, hopefully, give a good order of magnitude estimate in higher dimensional spaces, as well. 



STABLE BLACK HOLES AS COMPONENTS OF DARK MATTER 



If black holes become stable after approaching a nonzero minimal mass then their number can only be reduced 
by collisions with each other, when the mass of one of the black holes is converted into radiation. Therefore, most 
primordial black holes and those produced later by cosmic rays should still be around at the present time. For large 
black holes the higher dimensional Lovelock terms become irrelevant and the Schwarzschild approximation is valid. 
Thus, the spectrum of black holes produced after inflation is unaffected by them. Hawking and Page [57] have shown 
that primordial black holes lighter than 10^^ g would have decayed by now. However, in odd dimensional Lovelock 
theory they would have all ended up in the near-minimal mass state. These black holes, along with small black holes 
produced later, would form part of dark matter. In the current section we investigate the effect of such a component 
of dark matter and derive bounds on their minimal mass, Mmin, and their relative abundance compared to other 
components of dark matter. 

First we show that for a very wide range of the Planck mass, Mp, and of minimal mass, Mmm, the passage of 
non-relativistic black holes through matter follows a simple pattern: A minimal mass black hole accretes a parton or 
an electron increasing the mass and the Hawking temperature of the black hole. The accretion is followed by the rapid 
decay of the black hole towards the minimal mass state, very probably before the next accretion takes place. This 
sequence of events allows us to derive a recursion relation for the velocity distribution and an integral equation for the 
limiting distribution. Not surprisingly, in view of the central limit theorem, the solution of the integral equation is a 
Maxwell distribution. It is interesting that the average velocity in the limiting distribution, depending on Afmin, can 
be higher than average velocity of incoming dark matter particles. Using the velocity distribution we show that in 
neutron stars the near-stable black holes are confined to a sphere of radius r < 100 m where they reach an equilibrium 
density in a few decades after the formation of the neutron star. In equilibrium, the number of black holes annihilated 
in pairwise collisions is equal to the number of black holes captured by the neutron star from dark matter. We will 
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also show that black holes cannot be captured by ordinary stars, if Mmin = 0(1 TeV), as the velocity of most of black 
holes exceeds the escape velocity. However, if Mmin, Mp > 100 TeV, then their velocity is smaller and the captured 
black holes accumulate inside the sun, as well. 

At the end of this section we will investigate the experimental signature of stable black holes passing through 
experimental devices. The signature is quite universal. It is in large measure independent of the minimal mass, up to 
Mmin — 10^ TeV, and of the velocity of the black holes. 



Decay of near minimal mass black holes in vacuum 

When we calculate the decay of black holes we can use the Schwarzschild metric for black holes with M » M„ 



while (11) can be applied to black holes with M < 2Af,„in. Considering that standard model particles live on a 3 



brane the radiation happens mainly on the brane 28 . Then the decay rate of near minimal black holes in vacuum is 
where g is the number of effective degrees of freedom. In 7 dimensions an extra factor of 2/9 should be included on 



the right hand side of (12 1. Note that, unlike for Schwarzschild black holes, in a theory with stable black holes the 



contribution of emission into the bulk is negligible when M is close to Mmin- The contribution of bulk radiation to 



( 12 1 would have higher powers of rh and , both vanishing when rh — >■ 0. 



Integrating (12) gives 



M(t) - Af„,i„ + , t^^, = , (13) 

where AM = M{Q) — Af,„in is the energy gained by the absorption of a parton. The half-time (time required to reach 
M(ii/2) = (Af(0)+M^in)/2 is 

ii/2 ~ 0.06^^ f • (14) 

When non-relativistic black holes travel through matter they absorb partons from atomic nuclei. Since they cannot 
leave colored matter behind, they must pull out an antiparton from the chromosphere of the black hole, as well. 
Altogether, they only gain about a mass of AAf = O(0.5GeV). Note that AM is, in a very good approximation, 
independent of Afmin. When, following the accretion process, owing to the limited available total mass, AAf, they 
can only emit light degrees of freedom, mostly photons, electrons, positrons, neutrinos and antineutrinos, with an 
occasional muon or pion. Then the effective number of degrees of freedom, including greybody factors j29) is 5 ~ 



23 2/3. Note that for 7 dimensional black holes the prefactor in (14) should be decreased to 0.03 



Applied to 5 dimensional black holes we obtain the following expression for 



^^^^ 0.013^ (^) .5.2xlO^TeV-^ (15) 



Note that here and all in subsequent equations Mp and M-^^-^ are measured in units of TeV and are therefore 
dimensionless. 



As we will see below, the decay time (15) is much shorter than the average time between collisions. Then it follows 
that the passage through ordinary matter is a series of accretions of partons each followed by almost instantaneous 
decay back to near minimal mass. 



Cross section of TeV scale black holes 



For the sake of definiteness we will consider black holes in a RS I scenario. A classical calculation of the Bondi- 
Hoyle-Littleton accretion can be performed by considering the geodesies of pointlike incoming particles at impact 
parameter h. If the impact parameter is smaller than the compactification radius then the absorbtion cross section is 
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defined by the maximal Trb"^^^, where bmax is the maximal impact parameter at which the projectile is captured by the 
black hole. For this calculation we must use the D dimensional Einstein-Gauss-Bonnet metric in the neighborhood 
of the black hole. We obtain for Z? = 5 

a = Aa^.v-Hlf (2 + ^) = g^^^minMp^ (2 + ^) ' 

where v is the relative velocity and 7 is the corresponding Lorentz factor. Since we are interested in order of magnitude 
estimates we use classical cross sections throughout this paper. 

The increase of a with decreasing v is due to the increase of the maximal impact parameter at which the particle 
is captured. However, at large impact parameters the 5 dimensional metric is no longer applicable, because at large 
distances from the black hole the metric becomes effectively 4 dimensional, with standard 4 dimensional gravitation 
strength, which hardly causes any deflection of particles. The transition between the 5 and 4 dimensional behaviors 
happens at r ~ where rrij, is the radion mass. The largest lower bound on the radion mass, 120 GeV, has been 
derived in [301 ■ We are not aware of upper bounds. 

To incorporate the low velocity cutoff imposed by the radion mass we modified our calculation of the cross section by 
using the simple minded assumption that the metric reduces to Minkowski metric at distances larger than Tmax = iti~^. 
Then at nir << Mp and v << 1 



(17) 



For larger black holes the cross section can be calculated from the Schwarzschild metric. One obtains [3T] 

a = Anfil, (18) 

where /iu is given in Q. In particular, for D ~ 5 

a = '^MMp^ (19) 



Mean free path of black holes 

The mean free path in a medium is given by / — (crp)^^, where p is the number density and a is the collision cross 
section. To get a rough estimate of the mean free path we will use the absorption cross section. We will discuss 
neutron stars first. 

Nuclear matter has the parton number density, with about 3 partons per neutron 



p = 4 X lO^^^TeV^ 



Then using (16) and (17), we obtain 



1.1 X lO^TeV"^ < L, 



— < 3.8 X 10" TeV-i 
pa 



M„ 



(20) 



(21) 



where the upper limit is given by the minimum of ( 16 ), at 7 — > 00 and the lower limit by (17) at the minimum radion 
mass [30] =120 GeV. However, if Mp > 10 TeV, then the velocity dependent classical formula gives a larger lower 
bound then the one obtained using the minimum radion mass ((w^) = x AM will be calculated in the next 

subsection) : 



3.8 X 10^' TeV" 



I 



= 2 X 10« TeV" 



(22) 



The lower bound (22) ensures that black holes undergo many collisions in nuclear matter (or in stars), except in the 
case when our mass parameters, Mp and Afmin are at the upper end of the range of interest, near Mp ^ Mmin > 10^^. 
Then the mean free path becomes larger than a few km and most black holes undergo a small number of collisions 
only. Since the number of collisions inside the star largely depend on the total mass only, a similar statement applies 
to a star, like our sun, as well. 
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The corresponding bounds on the mean free path in water are 



Af3 (2.9 km 

600km— >C > <^ Ml (23) 

JUmin 0.3 km X ■ 



where the first lower bound is obtained using cross section (f7), while the second using (16) at low v. These numbers 



will be important when we will discuss experimental constraints on Afnjin and 5, the fraction of dark matter in the 
form of microscopic stable black holes. 

While passing through earth black holes undergo a large number of collisions, unless Mp > is very large. The limits 



for water, (23 1, should be divided by 5.5. Even though they collide with nuclei many times inside the earth, the earth 
cannot capture black holes, because, as we show in the next subsection, their average velocity, {v) = O(10~^c) is 
orders of magnitude higher than the escape velocity. 

Let us compare now the decay time of black holes with the collision time in nuclear matter. The average time 
between collisions is tc = (n/v. In the next subsection we will determine the average speed of black holes, dependent 
only on having a large number of collisions and largely independent of the matter density. We will obtain v ~ 
0.5 -y/ AM /Minin — 0.01 Mj^jj/'^. Using (21) and taking the accreted energy AM ~ 0.5 GeV in each colhsion, we 



obtain the following expression for the ratio of collision time and decay time in nuclear matter 

m6 / fl.SxlO^^ 
6.8xl0^^>-^> (24) 

Taking the root of the product of the two lower limit, setting Mp ~ Mmin, and using the very conservative value 
X = 0.2 we obtain a lower limit that is independent of Mp, namely > 9 x 10^. Even in the center of the 

neutron star where the density is 2-3 times higher than that of nuclear matter the ratio is much larger than 1. 



The important message in ( 24 1 is that the passage of black holes through any medium, including the core of neutron 
stars, follows a simple pattern. After a M ~ Mmin black hole collides with and absorbs a parton its excess mass 
evaporates long before the next collision to become a M ~ Mmin black hole again. This sequence will be important 
in the next subsection in which we investigate the velocity distribution of non-relativistic black holes close to their 
minimum mass in matter. 



Diffusion of black holes in matter 



We consider stable black holes arriving, as part of the dark matter, at a star or a planet, in which they undergo 
multiple collisions. For simplicity, we will consider the star as infinite matter. If the temperature of the star is low 
compared to the kinetic energy of incoming stable black holes one may naively think that the black holes lose velocity, 
until their kinetic energy reaches the ambient temperature. In fact, this is not the case at all because the inelastic 
collision and subsequent decay of the black holes may result in, depending on the mass of the black hole, an increase of 
the average velocity. The black holes are not in equilibrium with the surrounding matter, because of the long collision 



time compared to the decay time, as shown by (24). When the black hole accretes a parton then, as an average, it 
slows down slightly in the star rest system. However, after accretion it decays soon approaching its minimal mass 
state by emitting a small number of particles and gaining extra velocity from the recoil. After multiple collisions 
the addition of these recoil velocities results in a random velocity distribution, which is independent of the original 
velocity distribution of incoming dark matter particles. The average velocity depends on A/min but as we will see 
below if Minin ^ 1000 TeV then it is higher than the average velocity of dark matter particles. 

Let us consider the first step of the repeating two-step process, the accretion. As black holes accrete partons of 
energy AM then, as an average, their momentum is conserved in the star frame but they lose approximately the 
fraction AM/Mmin of their velocity. However, in the second step, when they subsequently radiate essentially all of 
the accreted AM energy they also gain momentum 

n 

P+ = -E*' (25) 



where qi are the momenta of emitted particles. The n particles are emitted isotropically. Since the Hawking temper- 
ature Th » AM the individual particle momentum distribution can be approximated by phase space. Neglecting 
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the mass of emitted particles the momentum gained by the black hole after the end of Hawking radiation is 

p / " \ / n \ n 

f(-\p)=N-' / ,5 p + ^ql h AM-^q, md^ft, (26) 



where qi = \qi\ and A'' is a normalization constant. Note that n, the number of particles emitted, is small, (n) 2, 
as it can be ascertained using the infinite temperature limit of the Bose-Einstein and Fermi-Dirac distributions. In 
particular. 



1 

N 



(27) 



where N — AAf'^7r[log(4096) — 17]/36. The distribution f^'^Hp) is slightly more complicated, containing polyloga- 
rithms. 

Now the question is whether after successive accretions and subsequent Hawking radiations the momentum dis- 
tribution approaches a limiting distribution or not. Originally the speed of the black hole is about v ~ lO^^c, a 
typical speed for dark matter particle. The change of the speed of the black hole due to the emission of particles is 
5v = 0(p/M„iin) = 0{AM/Mtnin)- This is of the same order of magnitude as the original velocity, at least if the mass 
of the black hole is M = O(lTeV). Even if the black hole has mass Mmin >> ITeV it loses memory of the original 
direction, as it undergoes a very large number of collisions. Thus, its distribution becomes isotropic (discounting 
the effect of the gravitational force of the star). After a few collisions the speed distribution of the black hole in 
homogeneous infinite matter approaches an asymptotic form determined by the integral equation 

9iv) - ^ / 5 (^^[1 + AAf/Af,„i„] - g7Af,ni„) f{q)d^q (28) 

where the multiplier 1 + m/Afniin is the factor by which an accretion slows down the black hole and f{q) is the 
normalized weighted average of distributions f^'^^q) over k. 



As expected from the central limit theorem, the solution of (28), in leading order of AAf/M,„in, is a Gaussian 
distribution, 

/3AMM^i„y/' / 3AMM^int;2\ 
2n{q^) ) m j' 

where 

{q^) = J q'f{q)d^q = xAM', (30) 

with a; < 1. For massless particles the value of Xn {x for n decay products) is a;„ = (1 + 3(n — l)/4)~^. 
The expectation value of the squared velocity is 

(^')=^^^^'<^^~'^M-l (31) 

-'''-'min 

One would expect that the number emitted particles is small, so 1 > x > 0.2, leading, for My^in — 1 TeV, to an 
average velocity, Uav ^ 10~^c. In other words, the multiple collision process accelerates the minimal mass black hole, 
provided Af^in < 1000 TeV. 



Eq. ( 29 ) is a Maxwell distribution with kT < 200 McV. It is interesting to observe that the temperature of the gas 
of black holes is independent of the temperature of the medium. The two different temperatures are maintained by 
the peculiar nature of the interaction of black holes with particles in the medium. 



Near minimum mass black holes in stars 



The calculation of the average velocity in the previous section has assumed that no external forces act on the black 
hole. This is not true in stars and particularly in a neutron star in which the gravitational field is enormous. Using a 
radius of R =10 km, and mass of A/„ = 1.5 x Afgun we obtain the gravitational potntial, g, as 

ff=^=Q.2»^<2.5xlO-^ (32) 

XL 2 
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The velocity of black holes is then much smaller than the escape velocity and the black holes start to move towards 
the center of the neutron star. At the center, however, the gravitational field vanishes, so the black holes rapidly reach 
the velocity distribution described in the previous subsection. They will be confined to a sphere the approximate 
radius of which is obtained from equating the average random velocity with the escape velocity on the surface of this 
sphere. To determine this radius we need a model for the density distribution. We choose a normalized quadratic 
density distribution of 

Other choices would slightly change the final result. Then near the center the gravitational potential is 

10 

g(r)^M„G-^. (34) 



Then using (31 1 and (32) we obtain the the radius, Tc, of the aforementioned sphere from 

i,2\ ^2 



We obtain 



^=^^^ = 2]^ (35) 



-llOm^Mjf. (36) 



Inside the sphere of radius black holes annihilate by colliding with each other. In an equilibrium state the number 
of annihilated black holes must equal the number of black holes captured by the neutron star from dark matter and 
from those produced by cosmic rays at the surface of the star. 

If minimum mass black holes constitute a fraction S of dark matter, then using the dark matter density of 300 GeV 
m~^, a radius of 10 km , and a mass of 1.5 times that of the sun we obtain the following capture rate of black holes 
traveling at ~ 270 km/s 

iVeap = 2.5 X lO^^s-i SM;^!^ (37) 
The number of black holes undergoing two-body collisions in the central sphere per unit time and unit volume is 

Na = pIh^^v (38) 



Then using the Maxwell distribution ( 29 1 the average relative velocity of two black holes is 

aT/ AM 
37r \Mmii, 

Now after attaining equilibrium VNa = Neap so the density of black holes in the central sphere is 

= ISm-^ = 

vVav x^/a/ 



(39) 



PBH = = 18 m-^ -^^=^M^J,l (40) 



where V is the volume of the central sphere confining the black holes. The cross section is measured in m^. 

Now we can use our three bounds on a to find bounds on the density of black holes in the central region of a 
neutron star. Using (16) in the extreme relativistic limit we obtain an upper limit, while using (16) with (w^) taken 
from (31) provides one of the lower bounds, and finally (17) provides another lower bound, as follows 

1.8 X lO-m-^M^ > PBH > I' ^ 10"--^\/i<'^^-n^ (41) 

^2.9x 10i«m-3My' 



min ■ 



These limits can be compared with the density of partons in the central region, pp ~ lO'^^m"'^. 

How long a time is needed to fill a newly created neutron star with black holes at equilibrium density? The total 
number of black holes in the central region divided by the rate of black hole capture provides two lower bounds on 
ifiii, obtained from (41 1 

PBH ^ J il)^^^ years 



- A^cap- l3^A4/^M-i^,years 
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It is worth noting that if Mp and Mmin are close to M^'^ ~ 2.4 x 10^^ TeV then the process of "filhng" the neutron 



star with black holes, which, according to (36), are now concentrated at the very center of the neutron star, takes tens 
of millions of years. This would imply a change in the surface temperature on the scale of tens of millions of years, as 
well. An investigation of time dependence of the surface temperature and its possible observation in neutron stars has 
been performed assuming that WIMPs constitute dark matter [TU]. Black holes as components of dark matter, could 
effect the surface temperature of neutron stars considerably less than WIMPs. The mass of the black hole produced 
from the collision of two black holes may be substantially less than the sum of the masses of the colliding black holes 
j20j . The lowest possible resultant mass is obtained from adding entropies, rather than masses. In the case of near 
minimal mass stable black holes that would imply a third black hole of near minimal mass, as well. The rest of the 
energy, A^E ~ Mmin, would be dissipated in the form of gravitational radiation, which would escape the neutron star. 
In such an extreme case the annihilation of stable black holes would not lead to the heating of stars, at all. However, 
if the collisions are not entropy conserving then the masses of the black holes may simply add. In the latter case a 
large fraction of the mass would be radiated in the form of elementary particles. We intend to return to this question 
in a future publication. 

We will consider ordinary stars briefly. In a star like our sun, or in any star, for that matter, black holes undergo 



a large number of collisions to reach the equilibrium velocity distribution ( 29 1 . Owing to the large size of the sun or 
other ordinary stars, the surface gravity is smaller then the average squared velocity, provided the minimum mass is 
smaller than Mmin < 100 TeV, 



MsunG ^ „ ^„ii, , n9 1 



g = = 1.9 X 10"(m/s)2 « -(u^) ~ 1.8 x 10"a; M^ijm/s)^. (43) 



However, if M^in > 100 then black holes get captured in the sun, as well. In fact, the number of black holes collected 
in a central sphere is much larger than those in a neutron star, owing to the sun's much larger cross section to capture 
dark matter. We intend to address this question along with the question of dissipation of heat produced from the 
annihilation of black holes from stars |19j in a future publication. 



Signature of stable black holes traversing matter and bounds on the minimal mass and abundance 

As we pointed out in the previous subsection, stable black holes contributing to dark matter travel through visible 
matter in an unusual manner. Suppose at a given instant the black hole is close to its minimal mass with nearly zero 
Hawking temperature. Then colliding with an atomic nucleus, or an electron it can absorb a parton or an electron 
and increase its mass and Hawking temperature. Note that stable black holes near the minimum mass have positive 
specific heats. The mass of the black hole increases by the energy of the accreted object, AM << Mmin, while, as 



shown by (11 1, its temperature increases to Th ~ V Mmin AM. This implies that the increase in temperature is larger 
by at least two orders of magnitude than the increase of energy. The black hole becomes a very hot object with very 
low heat capacity. 

It is difficult to calculate the average particle number and momentum distribution in this strange system. However, 
owing to the fact that the total energy to be radiated isotropically is in the range of AAf ~ 0.1 — 0.5 GeV the 
decay products are restricted to light particles, mostly photons and electrons and neutrinos, with a possible muon or 
pion. A certain amount of gravitons are also emitted, further decreasing the missing energy [32] . Due to phase space 
considerations lower multiplicity is probably preferred, so most of the events should contain 1-3 visible particles. As 
the recoiling black hole takes up momentum, the total momentum is not balanced even if we disregard the momentum 
taken away by gravitons and neutrinos. In many ways the events should be similar to the decay of a proton-like 
object with a mass of less than a GeV. Therefore, the ideal experiment searching for stable black holes is in large 
underground detectors. 

As in proton decay experiments, the most important background signal is the inelastic interactions of atmospheric 
neutrinos. In fact, the frequency of i? < IGeV u^, De events serves a rough upper bound on the frequency of stable 
black hole collisions. In the Minnesota proton decay experiment '33] 69 neutrino events, resembling proton decay 
were found in 3300 m'^ fiducial volume in 80 days. As stable black hole emission events also resemble proton decay 
events, we may use the frequency of those neutrino events as an upper bound on black hole evaporation events. This 
amounts to an upper bound on the event rate of A^e = 3 x 10~^ m~'^s~^. Compare this with the limits on the rate of 
predicted stable black hole events, using our lower and two upper bounds on the cross section : 

, . 1 f6200(5Mpim-3s-i, , ^ 

162 5 Mp3 m-^s-i < ATbh < < % ^ i 44 

P _ - |3 7 X 105iMminMpV-3s-l ^ ^ 
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Then the lower bound on the black hole decay rate provides a lower bound on the Planck mass that must be satisfied 
so that we would not contradict the observed rate of neutrino events resesmbling proton decay 



5 Mp^ < 2. X 10-". (45) 

The upper bounds on the black hole decay rate provide two further conditions. If either of these are satisfied then 
constraint on the event rate is certainly satisfied. The two conditions are 

S M,:^i^ < 5 X IQ-^^ (46) 

and 



Then (45) implies that if the dark matter consists solely of black holes then Mp > 3700 TeV. The bounds (46) and 



(|47|) imply that the event rate is certainly low enough if Mmin > 2 x 10^^ -peV or M^^M^^H'^ > 1.1 x ^lO'^ TeV. 



If Mp ~ A/jjiin then the lower bound on Mmin is somewhere in the range 

3.7 X lO^TeV < M^in < 1-1 x lO^TeV. (48) 

Black holes with a minimal mass below the lower bound would surely have been detected. If the minimal mass is 
above the upper bound then the black hole would certainly have not been detected. 

Unless only a small fraction of dark matter is made up from stable black holes the minimum mass of black holes 
must be much larger than 1 TeV. If the mass of black holes is very large, Mmin > 10^° TeV, then the decay event is 



not any more localized. Using ( 13 ) in 5 dimensions we obtain for the expected remaining mass as a function of the 



distance, L, starting from the point of accretion 



AM(L) ~ 2.4 X lO-i^TgV (L + Lq)"^/^^^, (49) 

where L is measured in meters and 

Lo = 2.310-"m-^ (50) 

Fig.l shows the dependence of the average energy emitted by the black hole per meter on L at five different choices 
of M — M^^^/Mp. The initial excess energy has been chosen to be AM = 0.5 GeV. 

We can conclude from Fig.l. that if M = M^p'' = 2.4 x lO^^ TeV there is a significant chance of observing a series 
of emissions of 7 rays and energetic electrons at large distances from the point of accretion. At lower values of M, 
certainly for M < 10^° TeV most of the energy is deposited within a very short distance and only some hard X-ray 
photons can be emitted up to approximate distance of 10 m from the point of accretion. 

The spatial resolution of the Super-Kamiokande detector for sub-GeV particles is around 0.3 m [33]. We can see 
that at M < 10^ TeV the decay appears to be completely localized at the point of accretion. 

The Super-Kamiokande detector is eminently suitable for discovering stable microscopic black holes, provided the 
abundance of black holes in dark matter is large enough and Mp is small enough. The approximate limits on these 



quantities are shown in Fig. 2. along with the upper bound on the flux (46) and (47). A further limitation on the 
mass of stable black holes comes from accelerator experiments at Fermilab or the LHC. All these limits are plotted 
in Fig. 2. 

One may consider the scenario of an eleven dimensional theory, in which there is a lower bound on the black hole 
mass in the presence of the Lovelock term L5, compactified above the Planck scale oi Mp < 10^*^ TeV. Then, if the 
minimal mass is of the same order of magnitude, all the bounds are satisfied, even if all of dark matter is constituted 
from black holes. Unfortunately, using (44), if Afmiu ^ Mp > 10^ TeV the event rate at Super-Kamiokande drops 



below 1/year and observing an accretion and decay of a black hole becomes not feasible. One may significantly 
increase this upper limit if one also searches for series of 7 emissions that a black hole emits after it accreted a parton 
in the rock outside the fiducial volume of the detector, as shown by Fig. 1. 
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FIG. 1: The average deposited energy per meter for M = M^p ^ (solid line) , PI = 10^^ TeV (dashed), M = 10^" TeV (dotted), 
M = 10** TeV (dash-dotted), M = 10^ TeV (bottom solid hne) 
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FIG. 2: The 5-M plane: (1) The region left of the dotted line is excluded on the basis of accelerator experiments, (2) The 
region above the dash-dotted line is excluded by hmits on the flux in proton decay experiments, (3) The region between the 
dash-dotted and dashed lines represents stable microscopic black holes that are amenable to discovery at Super-Kamiokande, 
(4) In the region under the dashed line the event rate is too low for discovery with currently available detectors, (5) the 5 > 1 
region above the solid line is excluded because of astronomical limits on dark matter. 



RELATIVISTIC BLACK HOLES IN A MEDIUM 



In low scale gravity, 2.4 x lO^^ TeV = APp' > Mp > 1 TeV, high energy particles can create black holes. If 
the primary energy is i? >> Mp/nip, such as may be in energetic cosmic rays, then the produced black holes have 
relativistic speed and they may or may not slow down and evaporate while passing through matter. The mass and 
speed of the black hole changes due to the absorption of matter particles and Hawking radiation. We will investigate 
below what are the conditions for the black hole, undergoing multiple collisions, to slow down and evaporate before 
accreting all matter available to it |35j . 

Since in a brane- world scenario all particles save the graviton and Kaluza-Klein modes are confined to the brane and 
the Kaluza-Klein modes contribute little [55], despite the fact that the black holes are higher dimensional, we write 
the equation for the time dependence of the mass in the rest system of the black hole using 4-dimensional quantities. 
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The change of mass of the black hole as described in the rest system of the black hole is 



dM 



2 " ^ 2ti4 



(51) 



where g = X^s-'^sSeffi where and g^g are the gray body factor and effective number of degrees of freedom for 
particle species s [37]. c is the effectiveness of turning the energy of a parton into the mass of the black hole. Stefan's 
constant per massless degree of freedom is 7r^/120, p is the mass density of the medium and 7 is the time dependent 
Lorentz factor, r/j is the radius of horizon. The multiplier 7^ should be included in the absorption term because on 
one hand the black hole sees a Lorentz contracted distribution and the other hand the energy of each particle of mass 
m absorbed is 7717. Another way of stating this is that the accretion rate in the star's rest system is 



dM 



= ca p Aw, 



(52) 



where Aw is the relative speed. Now transforming into the system of the black hole while tp — > j'^tp and the relative 
speed is unchanged. We defined 



where p, the radius of horizon, is related to the mass as in Q. Th is the Hawking temperature. 



Th 



2a 



27r(/z2 + 2a) ' 



(53) 



(54) 



Since 7 is also time dependent we need another equation to find both p and g. We can readily obtain such an 
equation if we consider the process in the rest system of the medium. Then using momentum conservation in the 
rest frame of the star and assuming that there is no momentum transfer to matter not absorbed by the black hole we 
obtain 



d{M-iv) 
Jt 



(55) 



where we assumed that the c.m.s of evaporating particles moves with the same speed as the black hole. 
Combining (51) and (55) we obtain an equation for the Lorentz factor 



(^7 
It 



2 _ -|N3/2 



M 



cpiY - 1) 



In the Schwarzschild limit {M >> a Alp, where Mp is the Planck mass) the ratio 



M 



(56) 



(57) 



is a constant and the equation for 7 decouples. Introducing the notation x = 4:cp/{3Mp) and using the initial 
condition, 7 = 70 at t = 0, the solution of (56 ) is 



1 + xtvo 



(58) 



p{y) can be obtained as a functional of 7 from (51) in the Schwarzschild approximation, with initial condition 
M'(0) = Mo = A^o/(37rM3), 



[h{z)]-^dz 



(59) 



where P = gir ^A/p'^/180, y — xt, and 



h{y) = exp 1^ [7(z)]2dz| = 7o\/(2/«o + 1)^ - "o- 



(60) 
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Calculating the integrals in (591 we get 



= 70 [{xtvQ + ly 



211/2 



4 



log 



(1 + Uo)(l +VoXt- Vq) 



%xvo \{l-vo){l + voxt + vo) 
Now there are two possible scenarios. The first one is when /i^ — > oo when i — > oo, if 



1/2 



Mo > 



log 



l + t-o 
1 - Wo 



/3 



lixvo 



log(47o')- 



(61) 



(62) 



In that case the mass of the black hole grows indefinitely, i.e it gobbles up all matter available. In the opposite case an 
Einstein black hole evaporates completely in finite time. The evaporation time can be obtained from finding the zero 
of the last factor of (61 1 . For a Lovelock black hole, which approaches a finite mass as t — )■ oo, we must slightly modify 



(611 at large times. The difference between Einstein and Lovelock black holes in no way modifies our conclusions 



concerning the fate of cosmic ray black holes in neutron stars. 

Using Mo7o = Eq, the primary energy, we obtain the following upper bound on the primary energy for the black 
hole to evaporate already in the crust of the neutron star (with density 10^ kg/m'^) 



Eo < 10' 



^TeVMl,. 



(63) 



According to current theory and experiment the cosmic ray spectrum cuts off at i?max = 5 x 10^ TeV due to the 
Greisen-Kuzmin-Zatsepin effect f36 . We estimate that a black hole created by a cosmic ray proton of E < i?max in 
the crust of a neutron star will lose its energy in approximately / ~ 1/a; ~ 1/im. 



CONCLUSIONS 



In more than 4 space-time dimensions Lovelock terms make a non-trivial contribution to the equations of motion. 
Putting it in stronger terms, there is no reason known to us why Lovelock terms should not appear in the gravitational 
action. Their presence significantly modifies the physics of small black holes. Some of these modifications, especially 
in odd dimensions, in which black holes have a lower bound on their mass, have been investigated in the past [T2] [TT] 
[T5] . At the lower bound odd dimensional Lovelock black holes, just like extremal black holes, have vanishing Hawking 
temperature [23] [52] [H]. At least semi-classically, at the end of their decay process, these black holes become stable 
in vacuum. 

If D is even, then after the compactification of an odd number of dimensions, we believe that small odd dimensional 
black holes can be caged without a substantial change of their metric and thus they retain their property of having 
a minimal mass and vanishing Hawking temperature. This could be ascertained in an expansion in the inverse of the 
compactification radius [25] [26] . Alternatively, in even dimensions, if the compactification radius of an odd number 
of dimensions is smaller than /i*^^^ black holes become black branes. Such black branes retain the combination of 
properties of having a minimal mass and vanishing Hawking temperature at the minimal mass. This will be the 
subject of a future publication. 

If black holes become stable upon reaching their minimal mass then their number density can be reduced only by 
collisions of black holes. Then by necessity, near minimum mass black holes form a part of dark matter. One major 
concern of this paper has been to find upper limits on the contribution of stable microscopic (of mass between 1 TeV 
and the four dimensional Planck mass) black holes to dark matter. We have provided an approximate upper bound 
on their abundance using existing data at underground experiments. In fact, we also show that using the extreme 
assumption that all of dark matter is constituted from stable microscopic black holes their mass is at least lO'' TeV, but 
probably more than 10^ TeV. In doing so we show that stable black holes traversing through underground chambers 
have a very peculiar characteristic signals, especially if their minimum mass is larger then 10^ TeV. In the range 
-^min < 10^ TeV their decay is local, often imitating a low energy {E < 0.5 GeV) neutrino interaction. However, at 
higher energies, especially at Mmin > 10^° TeV, after a localized partial evaporation they continue emitting lower and 
lower energy particles (mostly photons and possibly electrons) along a track as they progress through the detector. 
Unfortunately, if Mmin ^ Mp > 10^° TeV then the event rate at the Super-Kamiokande detector becomes so low that 
the observation the black hole evaporation process becomes very difficult. 

Black holes that form part of dark matter accumulate in a central region of neutron stars. There they annihilate 
by colliding with each other. Eventually they reach equilibrium with the incoming flux and then they produce heat 
though possibly at a lower rate than other components of dark matter 119 1 . In ordinary stars, however the surface 
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gravity is not sufficient to keep the black holes from escaping unless the M^iin >100 TeV, in which case the average 
velocity of black holes becomes smaller than the escape velocity. However, if the Mmin ~ Mp > 10^^ TeV the mean 
free path in black holes and in ordinary stars becomes comparable with the size of the star and most of the black 
holes are not captured. 

Among other aspects of the physics of black holes we investigated the accretion of relativistic (cosmic ray) black 
holes in matter. We find that all black holes created by cosmic ray protons below the Greisen-Kuzmin-Zatsepin limit 
[36j slow down and lose mass in stars, among others, in neutron stars. 

Finally, we return to a problem, which we already mentioned in the Introduction. Two further steps are required 
to complete our scenario, either of which could potentially destabilize minimum mass black holes. Suppose that D is 
the (odd) number of dimensions of Lovelock gravity, which we have considered in this paper. Since there are only 4 
non-compact dimensions in our world we still need to compactify D — 4 dimensions. 

The second problem concerns the existence of dimensions beyond D. The total number of dimensions of the world, 
2) > D. 2) may be even or odd. We must embed black holes into the J) dimensional space. In the next paragraph we 
will sketch a possible solution to these problems. 

Suppose the compactification scales oi T) — D and D — 4 dimensions are £ and L, where £ << i. Wc assume 
that D is odd and all Lovelock terms up to n = (Z) — l)/2 order appear in the action. The minimum mass should 
be small enough, such that black holes close to the minimum mass are caged in the D — A compact dimensions, but 
at the same time it should be large enough so that black holes are extended as black branes into the extra D — D 
dimensions. In other words, the condition for the existence of a stable minimum mass is that Ai(^) defined in Q is 
in between the two scales, i.e C£ £ < fi^^^ < c^L, where C£ and are 0(1) dimensionless constants. C£ is defined 
to be large enough to avoid the Gregory-Laflamme instability [5], while cl is chosen to be small enough to have the 
black hole safely caged in £> — 4 dimensions. The details of this scenario are left to a future publication. 
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